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Abstract 

The gravitational instability of a fully ionized gas is analyzed within the framework of linear 
irreversible thermodynamics. In particular, the presence of a heat flux corresponding to generalized 
thermodynamic forces is shown to affect the properties of the dispersion relation governing the 
stability of this kind of system in certain problems of interest. 
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I. INTRODUCTION 



The study of transport processes in plasmas leads directly to a set of partial differential 
equations that describes the evolution of the local thermodynamic variables relevant to 
its physical description. Linear irreversible thermodynamics predicts the coupling of all 
thermodynamic fluxes and forces involving all possible tensors of the same rank In this 
context, for multicomponent systems, the problem of gravitational instability in the presence 
of a heat flux associated to a density gradient (Dufour effect) is ought to be examined in 
order to address its implications regarding the Jeans instability criterion. The reciprocal 
effect, a mass flux due to a temperature gradient (Soret effect), does not affect the continuity 
equation for the total density of the system, and therefore does not alter the Jeans instability 
criterion. 

While lots of results regarding the effect of dissipation on the Jeans criterion have been 
obtained in several works to the authors' knowledge, this is the flrst time in which the 
Dufour effect is taken into account in calculations regarding gravitational stability. The 
numerical values for the transport coefficients used in this work have been recently obtained 
by the authors based on the Chapman-Enskog expansion used in order to solve Boltzmann's 
equation ^ Q|. Similar results have been obtained in terms of several formalisms, from the 
works done by Spitzer [sl, Braginski ^, Marshall I?] and Balescu 8|. 

The structure of this work is as follows. In section 2, the basic hydrodynamic equations 
are presented including the Dufour flux for a binary mixture of dilute gases. In section 3 the 
linearized system of equations for the fluctuations in the local density the local temper- 
ature 6T and the hydrodynamic velocity 6u are obtained together with the corresponding 
dispersion relation associated to it. Section 4 is devoted to the study of the explicit gravita- 
tional stability criterion for the binary plasma. The last section of this work is dedicated to 
a discussion of the results, addressing the conditions in which the various thermodynamic 
forces turn out to be relevant in the study of gravitational stability. 

II. BASIC EQUATIONS AND TRANSPORT COEFFICIENTS 

Consider a fully ionized system formed by species i = 1, 2, pi standing for the electron 
mass density and p2, corresponding for the ion mass density. The task to accomplish is to 
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determine the importance of the various dissipative effects on Jeans instability. Total mass 
conservation reads: 

^ + V-{pu) = (1) 
In Eq. ([T]) u corresponds to the hydrodynamic velocity of the mixture defined by the relation 

{pi + P2) u = pu = piui + P2U2. (2) 

The total density of the mixture is given by p = pi + p2- For a quasi-neutral plasma, 
the number densities are ni = n2 = |, where pi = mini and n = rii + n2, mi and 1712 
being the electron and proton masses, respectively. Also, Ui represents the average velocity 
corresponding to species i. 

n 

The balance equation for linear momentum reads [Jj: 

+ Vp + V ■ {puu + r^') = -pV^ (3) 

In Eq. (I3|), p represents the local pressure of the fluid , r*^' is the viscous contribution to the 
stress tensor and ip is the gravitational potential which, in Newtonian mechanics, satisfles 
the Poisson equation 

VV = -47rGp (4) 

Total energy conservation leads to a well known expression for the evolution of the local 

temperature of the system T, namely 

dT ( Be \ -> 
pc^— - j ■u = -V ■ (J[q]) - pc^,u -VT -pV -u (5) 

where pe is the internal local energy density. The heat flux vector J[q] in Eq. ([5]) contains a 
contribution due to a temperature gradient (Fourier law) and a contribution associated to 
a density gradient (Dufour effect): 

4] = -kVT - ^VT - ^Vp (6) 

where k stands for the heat conductivity and V is the Dufour coefficient. Equations ([3]) 
and (JHl), together with their corresponding constitutive equations relating thermodynamic 
fluxes and forces, form a complete set of non-linear partial differential equations which is in 
general difficult to work with. Nevertheless, stability analysis around constant equilibrium 
values may be readily performed by standard methods. This will be the subject of the next 
section. 
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III. LINEARIZED TRANSPORT EQUATIONS 



A first order perturbative stability analysis can be performed assuming that any local 
thermodynamic variable X in this system has a constant average value (X) and a fluctuation 
around it so that 

X={X)+ 6X (7) 

We can now rewrite the system given by Eqs. (JS]) and (JH]) in terms of the fluctuations. 
For simplicity, (u) = is assumed. The linearized continuity equation becomes: 

^ + W^» = (8) 

where 66 = V ■ 6u. Following the usual approach, the stress tensor in Eq. ([3]) is coupled 
linearly to the traceless symmetric part of the velocity gradient through shear viscosity. 
Thus, the linearized momentum balance equation becomes: 

(p) ^ + Ic'.W' i6p) + ip) (5T) - D^V' (56) = -AnG (p) 5p (9) 

where =^Vj V standing for the shear viscosity measured in international units {Pa — sec). 
In our dilute (ideal) plasma, bulk viscosity is neglected. We have also used the fact that, for 
an ideal gas, Vp = -C^Vp + VT where Cg is the adiabatic speed of sound and 7 is the 
heat capacities ratio (7 = 5/3). Therefore, the linearized energy balance equation is written 
as: 

m {p)c„ \ 2(T)I {pfc„2 >_ 3 ^ ' 

Notice that the Dufour coefficient V vanishes for a simple fiuid [Tl, l8||, but in principle affects 
the energy balance in the binary mixture. For our ideal gas, Cy = is the heat capacity 
measured in k = 1.38 x 10~^^^. n is measured in and V is given in — . The values 

KgK ' K Kms " j-ns i— i i— ■ 

of these coefficients for a dilute plasma have been obtained from plasma kinetic theory ^, l8||. 
They can be written as: 

.4Mm,2.oi.) (u) 

4 mi ^ 

^^5(n)(^ (12) 
2 nil 

Here r is a characteristic time given by 



where e is the electron charge, eo is he dielectric constant and ifj is the usual Coulomb log- 
arithm The shear viscosity coefficient rj is estimated by using the Eucken number for 
an ideal gas, 5/2 = n/rjCv We are now in position to address the importance of dissipation, 
including the Dufour effect, on the Jeans instability by means of the analysis of the corre- 
sponding dispersion relation. In order to derive the desired stability criterion, we perform 
a Laplace transform in time and a Fourier transform in space to the system (ISl fTOll . The 
corresponding equations become 

s6~p+{p)S6 = 6p{q,0) (13) 



AnG (p) - -Ctq') 6p + [D,q' + (p) s) 56 - -f^q'ST = (p) 56 (g, 0) (14) 



-^—qp + ^77^56 + 



If V \ , 



5f = 5f (g, 0) (15) 



2(p) 

Here, the symbol X = X {q, s) stands for the successive Laplace-Fourier transforms of the 
thermodynamic variable X. The wave vector is denoted by q and the corresponding Laplace 
frequency is s. 

It is interesting to notice that the first term in Eq. ([TSll has never been taken into account 
in the study of the Jeans instability. The term vanishes in a single component system and 
may become significant in the study of transport processes in plasmas, where the difference 
of masses betweenprotons and electrons is decisive while constructing the thermodynamic 
fluxes and forces [3|-[J|. 

The new dispersion relation reads: 

/ (s) = + as^ + /3s + 7 = (16) 



where the coefficients are: 



P = + , f;"" ,, q' + -^q' - (p) (18) 

V 2c„ (T) (p) c„ (p) J 

The absence of dissipative effects in Eq. ( fT6l ) leads directly to the ordinary Jeans wave 
number qj = J^^^§r^- 
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Figure 1: The function /(s) for values of q = qj — 6 for a number density n = lO^^m ^ and 
temperature T = lO'^K. 

IV. STABILITY ANALYSIS 

The behavior of the stability condition associated v^^ith the polynomial in Eq. ( flGl l de- 
pends significantly on its independent term vv^hich in turn contains dissipative eff'ects. We 
analyze two diff'erent cases to clarify this point, seeking critical solutions for wave numbers 
close to the ordinary Jeans wave number qj. 

First consider a system with n = lO^^m"^ and T = lO'^K. For these values of density 
and temperature, the cubic function / (s) in Eq. ( flGll . for a wavenumber q = qj — 6 (where 
6 is such that 6/qj <^ 1), has in general two critical points and the independent term is not 
negligible in the local scale here considered. In this case, the criterion for stability reduces to 
finding the threshold for which the maximum in the negative part of the s axis becomes zero. 
This can be seen from Fig. [T] which shows in dotted lines a case for which there is only one 
real root in the dispersion relation and one for which three real roots exist. The threshold 
that separates the ranges of q for which one obtains either pure damped or exponentially 
growing modes is around 6 = 6.6 x 10^^ which is shown in the solid line. 

Analytically, the maximum for s < is given by 

'- = -3 r- 
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Figure 2: The function /(s) for values of q = qj — 6 for a number density n = lO^^m ^ and 
temperature T = lO^K. 

and, since by varying q the whole curve evolves with /(O) virtually unchanged, the stability 
threshold is determined by = 0. 

The second case we shall examine corresponds to lower values of density and temperature, 
i. e. n = 10^^m~^ and T = lO'^K. The behavior in this case is quite different, as shown in 
Fig. [21 For these parameters, both a and the independent term are indeed negligible and 
thus, the critical value that gives the wave numbers is indicated by the disappearance of the 
two maxima. That is, the criterion can be found by imposing that the derivative of f{s) has 
no real roots at all, i. e. 

13 >0 

which, as clearly seen in Eq. ffTSl) simply reduces to the ordinary Jeans criterion in absence 
of dissipation. 

V. FINAL REMARKS 

The problem of the gravitational instability including dissipative effects has been studied 
within the framework of linear irreversible thermodynamics, including cross effects, for a 
dilute binary plasma. The dissipative effects enter both the hydrodynamic equations and 
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the dispersion relation for wave-like solutions of the linearized system. 

The wave number for which perturbations start growing exponentially in time does not 
differ much form the standard Jeans wave number and thus, the response of the system to 
different fluctuations wavelengths is similar to the case with no dissipative effects present. 
However, the qualitative behavior of the dispersion relation in terms of q is more compli- 
cated and has to be explored with care for each system. Cross effects are present in any 
multicomponent mixture and are enhanced by magnetic fields which are observed in many 
astrophysical systems These effects, which are predicted by irreversible thermody- 

namics and kinetic theory, are in general not taken into account while studying gravitational 
stability and further analysis of them should be performed in the future. 
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